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Abstract 

Given a set P of n points in the plane, we show how to compute in 0{n\ogn) 
time a subgraph of their Delaunay triangulation that has maximum degree 7 and is a 
strong planar i-spanner of P with t = (1 + \/2)^ * S, where S is the spanning ratio of 
the Delaunay triangulation. Furthermore, given a Delaunay triangulation, we show a 
distributed algorithm that computes the same bounded degree planar spanner in 0{n) 
time. 

1 Introduction 

Given a weighted graph G - {V, E) and a real number t > 1, a t-spanner of G is a spanning 
subgraph G* with the property that for every edge {p, q} e G, there exists a path between 
p and q in G* whose weight is no more than t times the weight of the edge {p,q}. Thus, 
shortest-path distances in G* approximate shortest-path distances in the underlying graph 
G and the parameter t represents the approximation ratio. The smallest t, for which G* 
is a t-spanner of G, is known as the spanning ratio of the graph G* . 

Spanners have been studied in many different settings. The various settings depend 
on the type of underlying graph G, on the way weights are assigned to edges in G, on the 
specific value of the spanning ratio t, and on the function used to measure the weight of 
a shortest path. We concentrate on the setting where the underlying graph is geometric. 
In this context, a geometric graph is a weighted graph whose vertex set is a set of points 
in yi'^ and whose edge set consists of line segments connecting pairs of vertices. The edges 
are weighted by the Euclidean distance between their endpoints. 

There is a vast body of literature on different methods for constructing t-spanners with 
various properties in this geometric setting (see [T] for a comprehensive survey of the area) . 
Aside from trying to build a spanner that has a small spanning ratio, additional properties 
of the spanners are desirable, e.g., planarity and bounded degree. 
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In this paper we consider the following problem. Given a set of points in the plane, 
the goal is to compute a bounded degree planar spanner of this set of points. Obtaining 
such a property enables us to perform routing algorithms that are known for plane graphs, 
e.g., p] and |3|. Bose et al. pJ| were the first to show the existence of a plane t-spanner (for 
some constant t) whose maximum vertex degree is bounded by a constant. Subsequently, 
Li and Wang [6] reduced the degree bound to 23. In [4], Bose et al. improved the degree 
bound to 17. Currently, the degree bound stands at 14, as shown by Kanj and Perkovic [s]. 
Here, we show how to reduce the degree even further. 

More precisely, given a set P of n points in the plane and the Delaunay triangulation 
DT(P) of P, we show how to compute in 0(n log n) time a subgraph of DT(P) that has 
maximum degree 7 and is a strong i-spanner of P with t - {1 + \/2)'^ * 5, where S is the 
spanning ratio of the Delaunay triangulation. This is a significant improvement in the 
degree bound from 14 to 7. Furthermore, given a Delaunay triangulation and a clockwise 
order of the edges, we show a distributed algorithm that computes a planar spanner of 
degree 7 in 0(n) time. Notice that we do not require the edges to be sorted by their 
length. 

Another result shown in this paper is Corollary |4.7| which states the stretch factor for 
a special case of the Delaunay triangulation. We hope that this result will help shed some 
light on the real stretch factor of Delaunay triangulation. 

Notation 1.1. Given a set of points P in , let DT(P) denote the Delaunay triangulation 
of P. 

Notation 1.2. Let Sp - {go, ■■■ ,Qk} be the set of neighbors of p in DT(P), labeled in 
clockwise order. 

Notation 1.3. Let Sp^q-^q. - {qk ^ Sp ■ qk is after qi and before qj in the clockwise order }, 
for angle ^{qipqj) < vr. 

For simplicity of presentation, in the rest of the paper we assume that in Sp^q.^q. the 
index of i is smaller than the index of j in the clockwise order. 

Notation 1.4. Let wedge Wp^g.^q^ = {{p,qk} ■ Qk e Sp^q-^g.} (see, Figure^. 

Notation 1.5. Let Psp^q^^q,, denote the path in DT(P) from qi to q^ restricted to points in 
^p,gi ,Qk ■ 

Notation 1.6. Let 5$^ denote the length of the path Ps^ ■ 

Notation 1.7. Let {p, (7min} be the shortest edges in DT(P) that is incident to p (qmin = 
minq^Sp\{p,q]\)- 

Notation 1.8. For a graph G - {V^E) and two points p,q ^ V, let Scip^q) denote the 
length of the shortest path between p and q in G. 
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Notation 1.9. For eachp^P, let Cp - {Cp^^ , Cp^ , . . . , Cp^^} denote a set of 8 closed cones 
labeled in clockwise order, with apex p and j angle, such that the boundary of cones Cp^ 
and Cpg contains the edge {p,qmin}- 

Notation 1.10. Let Dp^a,z denote the disk having p, a and z on its boundary. 




Figure 1: Wedge Wp^q-^q. and its edges depicted in solid lines 

The rest of this paper is organized as follows. In Section [2] we describe an 0(n log n) 
algorithm (Algorithm [T]) that bounds the degree of each point p e P by 7. In Section [s] we 
show that this algorithm gives a sub-graph of the Delaunay triangulation with bounded 
degree 7. Then, in Section |4] we prove that the stretch factor of the sub-graph resulting 
from the algorithm is a small constant. Finally, in Section [5] we show that given a Delaunay 
triangulation one can compute a distributed planar spanner of degree 7 in linear running 
time. 



2 Algorithm for bounded degree 7 

In this section we describe an algorithm that computes a bounded degree planar spanner. 
The approach we take to build such a spanner is to start with the Delaunay triangulation 
and then prune its edges to achieve the degree bound of 7 while maintaining a constant 
spanning ratio. To achieve this we ensure that for every edge of the Delaunay triangulation 
where we do not add to our resulting spanner, there is a spanning path in the resulting 
subgraph that approximates this edge. 

The algorithm consists of two main subroutines, BoundSpanner() and Wedge(). We 
start with the BoundSpanner() algorithm. The first subroutine (Algorithm [T]) is the main 
component of the algorithm. In the beginning of this subroutine, we compute the Delaunay 
triangulation (DT(P)) of the set of points and sort them in nondecreasing length order. 
We find for each point p its nearest point qmin in P] since {p, qmin} is an edge of DT(P), we 
use the DT(P) edges to find this edge for each point. Then, for each point p the orientation 
of the closed cones of Cp is defined according to this shortest edge {p,qmin}, such that this 
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edge is shared by cones Cp^ and Cpg. Next, an edge {p, g} is added to G if both p and 
q agree on it. A point p agrees on an edge {p, q} if every cone Cp. e Cp containing {p, q} 
is empty (i.e., Cp. n = for every cone Cp- e Cp containing the edge {p,q})- Notice 
that an edge could be contained in at most two closed cones, and exactly, two cones if this 
edge is on the common boundary of two cones. After adding an edge {p, we call the 
second subroutine {Wedge()) twice, once for point p and then for point q. In the second 
subroutine (Algorithm [2]) we add more edges to E that do not affect the degree of the 
spanner; however, they help us to bound the stretch factor of the resulting graph. 

Algorithm 1 BoundSpanner(P) 
Input: A set P of points in the plane 

Output: A planar t-spanner C = {P^E) with maximum degree 7 
1: Compute DT(P) ^ {P,Edt) 

2: Let L be a sorted list of the edges of DT by nondecreasing length 
■i: E ^ 
4: E* ^ Z 

5: Initialize Cp for each p e P /* with respect to edge {p,qmin} */ 

6: for each edge {p, g} e L (* in the sorted order *) do 

7: if (VCpj contains {p, q}, Cp- n E - 0) and (VCg^ contains {p, q}, Cq- n E - 0) 

/* Note: every edge can be contained in at most two adjacent closed cones */ 
then 

8: E ^ Eu{{p,q}} 

9: Wedge(p, q) /* calling subroutine Wedge() to check if some 

10: Wedge(g,p) edges (E*) are needed to be added to E */ 

11: E^EuE* 



Algorithm 2 Wedge(p, 

Input: Two points p and qi such that the edge {p,qi} e DT(P). 
Output: A set of edges E* to be added to the spanner C = iP,E) 

1: for every Cp^ contains {p,qi} do 

2: Let {p,qj} and {p,qk} be the first and the last edges in cone Cp^ (clockwise). 

3: E* <- E* u {{qm,Qm+i}} for each j < m < i - 1 

4: E* *- E* u {{qm,qm+i}} for each i < m < k - 1 

5: if (edge {p,qi+i} e CpJ and {qi+i t qu) and (angle ^(pqiqi+i) > 7r/2) then 

6: E* ^ E*u{{qi,qi^^}} 

7: if (edge {p,qi-i} e CpJ and {qi-i * qj) and (angle ^{pqiq-i-i) > it/2) then 

8: E* ^ E*u{{qi,qi_^}} 



Remark: Note that the output t-spanner C of P obtained by Algorithm [T] is a subgraph 
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of DT(P), therefore, it is planar and has a Unear number of edges. 



3 Bounded degree 

In this section we show that the degree of a point in the resulting graph is at most 7. We 



start by making two basic observations, then conclude with Lemma 3.3 



Observation 3.1. An edge {p,q} e Cp. is added to E during the loop in Algorithm^ step^ 
only if Cp. n E = 0. 

Observation 3.2. The first edge incident to a point p added to E during Algorithm^ is 
{p,Qmin}, thus, it belongs to two closed cones Cpi and Cp^ in Cp . 

Proof. Consider the time when edge {p, qmin} is considered to be added to E. Then, since 
it is the shortest edge incident to p in DT(P) it implies that all the cones in Cp are empty; 
thus, p "agrees" on adding edge {p,qmin} to E. Let Cq. be the cone in Cq that contains 
the edge {p, gmm}- Since the disk centered at p with radius \p,qmin\ is empty of points, it 
implies that cone Cq^ that contains {p,qmin} is empty. Thus, q "agrees" on adding edge 
{Pt qmin} as well, and the edge is added to E. m 

Lemma 3.3. The degree of spanner G constructed by the above algorithm is bounded by 7. 

Proof. Eight closed cones Cp are defined for each point p € P during Algorithm [Tj By 



Observation 3.2, there are two cones Cp^ and Cpg in Cp sharing a common edge. Consider 
the edges E^ incident to p that are added to E during Algorithm [l] (not including the edges 
added during Algorithm [2]) . Then each edge e e Ep is added to E only if the cone in Cp 
containing e is empty. Moreover, the first edge in Ep added to E shares two cones, thus 
l^^l < 7. Next, we show that the edges added during Algorithm [2] can be charged uniquely 
to empty cones, thus not increasing the degree bound of 7. Let {p, q} be an edge added 
to E during Algorithm [2] thus, there exists a point z such that the edge {p, q} has been 
added to E during the call Wedge(z,r). Moreover, this edge has been added in steps [3],|4] 
or during steps [6] , [8| 

• Case 1: The edge has been added during step|3]or|4j 

Let {p, s} be the consecutive edge to {p, z}, such that qi^ s. Since the edge {p,q} has 
been added to E during the call Wedge (z,r) it follows that the edge {z,s} is in the 
same cone (of Cz) as edges {z,p} and {z,q}. Thus, the angle ^{szq) < 7r/4 and by 
the empty cycle property of Delaunay triangulation angle ^{qps) > 37r/4. Therefore, 
there are at least two empty cones of Cp located between {p,q} and {p,s}. One of 
them is charged for the edge {p, q} and the second is left to be charged for the edge 
{p,s} if needed. 
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• Case 2: The edge has been added during step|6]or[8j 

In this case r - p. We know that the angle ^{qpz) > it/2, thus there is at least one 
empty cone c' of Cp located between {p,q} and {p,z}. Therefore, this empty cone c' 
is charged for the edge {p,q}. 

Therefore, the degree of every point p € P is bounded by 7. ■ 



4 Spanning ratio 

In this section we show that the spanning ratio of the resulting sub-graph is bounded. The 
empty circle property of Delaunay triangulations allows us to make two basic but crucial 
observations. 

Observation 4.1. From the empty cycle property of Delaunay triangulation it follows that 
each X e Ss,r,p is inside Dg.p^r- 

Observation 4.2. For qj,qi,qk ^ Sp^q.^q^ such that qi is between qj and q^ in the clockwise 
order, the angle ^qjqiqt > vr - ^qjpqt- 

Proof. Due to the empty cycle property of Delaunay triangulation, the point q^ lies inside 



the disk Dp^q-^q^ having p,qj,qk on its boundary (Observation 4.1). The angle ^qjqiqk is 
minimized when qi is on the boundary of Dp^q.^q^. In that case ^qjqiqk - n - ^qjpqk since 
the two angles lie on the same chord {qj,qk)- Therefore, ^qjq^qk > vr - ^qjpqk- ■ 

Observation 4.3. Let Dp^a,z be a disk having p, a, and z on its boundary and let (3 denote 
the angle -^(pza). Then, sm^/g) K^' ^11 length of the arc from p to a on the boundary 

of Dp^a,z (pa )■ 

Proof. Let o be the center of Dp^a,z and let r be the length of its radius, thus, angle 
^ (poa) - 2(3. By the law of sines, = 2r. Therefore, the length of the arc pa is 



2/3 2sin(^) sin(/3) ' 



Lemma 4.4. Consider a wedge Ws^r,p in DT(P) and assume that {s,r},{s,p} are the 
shortest edges in Ws,r,p incident to s (i.e., \sr\,\sp\ < \sx\ for all x e Ss,r,p\{r,p}). Then, 
^Ss,r,p{'^,p) ^ \rp\^i^, where a = ^{rsp). 

Proof. We prove the claim by induction on the rank of the angle a, i.e., the place of a in 
a nondecreasing order of the angles in DT(P). 

Base case: Angle a is the smallest angle in DT(P), thus, {r,p} e DT(P) and clearly 

^Ss.r,pif,p) < \rp\. 
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The induction hypothesis: Assume that for every a' < a the claim holds. 

The inductive step: If VFs^r,p\{{s, r}, = 0, then {r,p} e DT(P) and we are 

done. Otherwise, let {s,a} be the shortest edge in VFs^r.^p\{{s, r}, {s,p}}, i.e., {s,a} = 

^^^x^Ss,^A{vA^^'^\^ (s^^' Figured. 

Put a\ = ^{rsa) and 02 = ^psa). Since ai,a2 < ct, by the induction hypothesis, 
^Ss.r,ai'^,a) < \ra\^r^^, and 5s^,,,p(a,p) < lapj^jj^. Notice that S,^r,a £ S^.r.p and Ss,a. 
Sg f^p . Thus ; 



sin(ai) sin(a2) 



Note that by observation 4.1 a is located inside Dg^p^r- Let a' be the intersection point 
of -Ds,p,r and the extension of {s,a}. Since \sp\ < |sa|, ^{sap) < |, therefore, ^(a'ap) > | 
and |a'p| > \ap\. Symmetrically, we get \ra'\ > \ra\ and therefore, 

(?", p) < |ra I — - — - + \ap\- 



sin(ai) sin(a2) 



According to Observation 4.3, [ra^l ^-^"^^^ and lo'pljij^^ are the lengths of the arcs from 
r to a' {ra') and from a' to p (a'p) on the boundary of Dg^p^r, respectively. Moreover, 
I^Pl sin(a) length of the arc from r to p on the boundary of Dg^p^r^ which is the sum 

of ra' and a'p. Therefore, 



sin(ai) sin(a2) sin(Q;) 
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Claim 4.5. Let ^{pqr) he a right-angled triangle with hypotenuse {p,r). Then, 

TT TT 

{\pq\ + \qr\) < —\pr\. 



2^/2 

Proof. Let (3 denote the angle ^{prq). By the law of sines, \pq\ - |pr|sin(/3) and \qr\ = 
\pr\cos{/3). Therefore, 

vr vr 

(|pg| + |(7r|) = — ^|pr|(sin(/3) + cos(/3)) 



2V2 2V2 

<(-) 



-^^\pr\{^ + ^) 
2n/2 V2 V2 



2V2 

(*) The sum sin(/3) + cos(/3) is maximized when /3 = |. ■ 

Lemma 4.6. Let DT(P) be the Delaunay triangulation of the set of points P and let Ws,r,p 
he a wedge in DT(P), such that {s,r} is the shortest edge in Ws,r,p (\sr\ < \sx\ Vx e Ss,r,p)- 
Let r' be the projection of r on {s,p}. Then, 

TT 

where a = ^{rsp) < | (see, Figure^. 

Proof We prove the lemma by induction on the rank of the angle a. 

Base case: Angle a is the smallest angle in DT(P); therefore, {r,p} e DT(P) and clearly 

^Ss,r,p{'^,P) ^ \rp\- 

The induction hypothesis: Assume the claim holds for every angle a' < a. 
The inductive step: If Ss,r,p\{i^TP} - 0i then {r,p} e DT(P) and we are done. Otherwise, 
recall that from the empty cycle property of Delaunay triangulation it follows that each 
X e Ss,r,p is inside Dg^p^r- Let a e Ss,r,p be a point such that for every x e Ss,r,p\{i"}, 



\sa\ < \sx\. If |,sa| > \sp\ by Lemma 4.4 



and we are done. Otherwise, (|,sa| < \sp\), let a' be the projection of a on {s,p}. Denote 
ai - ^{asp) and 02 - ^{asr). Since ai < a, we can apply the induction hypothesis and 
get 

TT 

^Ss,rM^P) ^ ^^{\pO''\ + \a'a\). (1) 
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Moreover, by Lemma 4.4 



02 



sin(a2) 

a I 

sin(a) 



\ra\ 



2^2 



ra\ 



(since, a2 < a < —). 
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Therefore, 



'(2) 



TT 



2V2 



(1) _vr , , TT 



ra\ + — ^(Ipa'l + \a'a\) 



2V2 2V2 



TT 



\ra\ + \pa'\ + \a'a\) 



(2) 



The last inequality (*) is obtained by the following. Let b denote the projection of a on 
{r',r). Thus, \ba\ - \a'r'\ and \br'\ - \a'a\. Therefore, 



|ra| + Ipa'l + |a'a| 



\rb\ + \ba\ + \pa'\ + \aa'\ 
\rb\ + \a'r'\ + \pa'\ + \br'\ 
\pr'\ + \rr'\. 



Inequality (**) follows by triangle inequality, \ra\ < \rb\ + \ba\. 



/ a 




1 .■' 1^ 

■ '■ / 
'' '. / 


'a' 


\ 


r' 



Figure 3: Illustrating the proof of Lemma 4.6 



Claim 4.5 and Lemma 4.6 leads to the following corollary. 



Corollary 4.7. Let DT(P) he the Delaunay triangulation of set of points P and let s,r,p 
be points in P. If {s,r} is the shortest edge in Ws,r,p, then 

where a-^ (rsp) < j. 

Lemma 4.8. Let {s,r},{s,p} be two edges inCsjnDT(P) such that{s,r} has been chosen 
by Algorithm^ to be added to E. Then for every {s,x} €Ws^r,p,\sx\> min{\sr\,\sp\}. 

Proof. Assume, to the contrary, there is a point {s, x} e Ws^r,p such that \sx\ < min{\sr\, \sp\}. 
Let {s,w} be the shortest edge among all the edges in Ws,r,p- Since < \sr\, when {s,w} 
was examined by Algorithm [l| the cone Cg^ was empty from edges in E. Therefore, the 
only possible reason that could cause {s,w} not to be added to E is that E already con- 
tained an edge in the cone with apex w that {s,w} belongs to (w.l.o.g to Cwj)- Let {w,t} 
be an adjacent edge to {w,s} in C^j- Necessarily {t,s} e DT(P), and it is also in Ws,r,p- 
However, 

/^/^ TTVr 
^ {wtS) - -K - ^ [WSt) - ^ {swt) - ^ ~ 2 • — = — , 

we get \st\ < \sw\ in contradiction to the assumption that is the shortest edge among 

all the edges in Ws^r,p- ■ 

Claim 4.9. Let A(^rqp) be a triangle with ^{rqp) > then k\qp\+d\rq\ < k\rp\ fork > \/2d. 

Proof. Let q' be the point on {r,p}, such that \qp\ - \q'p\ (see, Figureji]). Since ^ (rqp) > 
then \rp\ > \qp\, therefore such a point exists. 

Since \rp\ - \rq'\ + \q'p\ - \rq'\ + \qp\, all there is left to prove is 

k\qp\ + d\rq\ < k{\rq'\ + \qp\), 

which is equivalent to 

d\rq\ < k\rq'\. 

Denote ^{pqq') - ^{pq'q) - /?, ^{qq'r) - 5, and ^(q'qr) - 7. Notice that /3 < |, therefore, 
I < (5 < vr, and f < 7 < f • 



By the law of sines, 



rq\ _ sin((5) ^ 1 ^ ^ 



\rq'\ sin(7)| sin(^) 
Thus, \rq\ < \rq'\\/2, which finishes the proof. ■ 

Claim 4.10. Let {s,r} and {s,p} be two edges in DT(P), such that \sr\ < \sp\ and the 
angle between {s,r} and {s,p} (^{rsp)) is less than 7r/4. Then, 

\sr\ + K(\rr'\ + \r'p\) < K\sp\ 

for K > , where r' is a point on {s,p} such that \sr'\ - \sr\. 

1 z sin(^ g ) 
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Figure 4: Illustrating the proof of Claim 4.9 



Proof. Put ^(r'sr) - a and ^(srr') - (3 - (see, Figure [s]); then by the law of 



smes, 



sin(o) 



Isr'l. Therefore, 



\sr\ + K(\rr'\ + \r'p\) < 



A*) 



sm(a) , , 
+ K{ . )' \sr\ + \rp\) 
sin(p) 

sin(a) /, I I I /|s 
+ K{ . )' \sr\ + \sp\-\sr |) 
sin(p) 

+ i^((^-i)M + M) 

sin(p) 

sin(^ 2 
sin(a) , 

cos(a/z) 
+ K(2sin(a/2) - l)|sr'| + \sp\) 
(1 + K(2sin(a/2) - 1)) + K\sp\) 
(1 + K(2sin(7r/8) - 1)) + K\sp\) 



K\sp\. 



The last inequality (*) follows from the fact that (1 + -R'(2 sin(Q/2) - 1)) is less than 

^OI K> i_2sin(f) - 



zero 




Figure 5: Illustrating the proof of Claim 4.10 



Lemma 4.11. The stretch factor of the resulting t-spanner of Algorithm 7] is (1 + \/2)^ • 6, 
where 6 is the stretch factor of Delaunay triangulation, i.e., t = (1 + \/2)-5. 
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Proof. Let G - {P,E) be the output graph of Algorithm [TJ To prove 
show that for every edge {s,p} e DT(P), (5g(s,p) < (1 + \/2) \sp\. We pro' 



the Lemma we 
prove the above by 

induction on the rank of the edge i.e., the place of the edge {s,p} in a nondecreasing 

length order of the edges in DT(P). 

Base case: Let {s,p} be the shortest edge in DT(P). Then, edge has been added 

to E during the first iteration of the loop in step[6| and therefore 6g{s,p) - \sp\. 
Induction hypothesis: Assume for every edge {r,q} e DT(P) shorter than {s,p}, the 
Lemma holds, i.e., 6G{r,q) < {1 + \/2)'^\rq\. 

The inductive step: If {s,p} e E, we are done. Otherwise, w.l.o.g. assume {s,p} e Cg. 
and {s,p} e Cpr, then, there exists either an edge {s,r} e Cs^ n E, such that \sr\ < \sp\, or 
an edge {p,r} e Cp^ n E, such that \pr\ < \sp\. Assume w.l.o.g. there exists an edge {s,r} e 



r\E, such that \sr\ < \sp\. By Lemma 4.8, for every x e Ss,r,p, \sx\ > min{\sr\, \sp\} - \sr\ 



Let {r,t} be the first edge in Pss,p,r7 ™d {q,p} the last. Note that all edges of Pss,t,q have 
been added to E during Algorithm [2] 

Claim 4.12. The edges {r,t} and {q,p} are shorter than {r,p}. 

Proof. If e € {{r,t}, {q,p}} is outside the triangle A(srp), by Observation |4.2| the angles 
'^(rtp), ^{rqp) > vr - ^[rsp) > and therefore, are shorter than {r,p}. Otherwise, it 
is bounded inside the triangle A(prr'), where r' is a point on {s,p}, such that \sr'\ - \sr\ 
(since \sr\ < \sp\, such a point exists). Therefore, this edge is shorter than {r,p}. ■ 

Applying the induction hypothesis on {r,t} and {q,p} results in: 

6G{r,t) < il + V2f\rt\ (1) 
SG{p,q) < {l + V2f\pq\- (2) 

Let {s,a} be the shortest edge in the wedge Ws,t,g (i-e., the closest point to s in 
5's,r,p\{r,p}). By Corollary 



4.7 



Ss.,.a ^ (3) 



2' 

TT 
2' 

Since \sr\ < \st\ and \sa\ < \st\, the angle ^{rta) facing towards s is less than vr. By 



^S,,,,, < ^\aq\. (4) 



Observation 
d = I gives us 



4.2 



we get that ^(^rta) > Applying Claim 4.9 on the triangle A(rta), with 



(l + ^2f\rt\ + ^\ta\ < {l + V2f\ra\. (5) 
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Therefore, 



< 
< 

<(1),(2) 
< 

<(3),(4) 



< 



(5) 



+ 5Gir,t) + 6Git,q) + 6G{q,p) 
+ {l + V2)H\rt\ + \pq\) + 6G{t,q) 
+ (l + x/2)2(H + H) + 55„,„ + 5s.,„., 
+ il + V2)\\rt\ + \pq\) + ^i\ta\ + \qa\) 

+ (l + v^)'(H + H) + |(H). 



There are two cases regarding the location of points q and a: 
• Case 1: Either point g or a is inside the triangle A(srp). 



\sr 




sr'\ 


sr 


< 


sq\, 



at s and with radius \sr\. Therefore, either point q or point a is located inside the 
triangle A(rr'p). 

1 



Since (1 + ^2)^ 



by Claim 



4.10 



we get. 



l-2sin(7r/8) ' 

\st\ + (1 + V2f{\rr'\ + \r'p\) < (1 + V2f\sp\. 
Therefore, it is enough to show that 



(1 + V2f{\ra\ + \pq\) + ^(|ag|) < (1 + V2f{\rr'\ + \r'p\). 

Observe the following two cases regarding the convexity of the polygon (raqp): 

— Case 1.1: The polygon (raqp) is convex. 
Since | < (1 + \/2)^, we get 

(l + v^)2(|ra| + M) + ^(H) < il + V2)\\ra\ + \qp\ + \aq\) 

< {l + V2)^{\rr'\ + \r'p\). 

The last inequality follows from the convexity of the polygon (raqp). 

— Case 1.2: The polygon (raqp) is not convex. 

The vertex that violates the convexity is either a or q, and the other vertex is 
inside triangle A(srp) (see, Figurejl]). Assume w.l.o.g. that a is the vertex that 
violates the convexity; then the angle ^ (raq) facing towards s is less than vr. 

(raq) > therefore, applying Claim 4.9 on the triangle 

(6) 



4.2 



By Observation 
A (rag) with d = f gives us 



(1 + \/2f\ra\ + ||ag| < (1 + y^Y' 



\rq\. 
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Thus, we get 



(l + x/2)2(H + |gp|) + |(H) <® (l + ^/2)2(H + M) 

{l + V2f{\rr'\ + \r'p\) 

The last inequahty (*) follows from the convexity of triangle ^{rqp). 





Figure 6: Illustrating the proof of Lemma 4.11, Case 1.2. In Figure (a) point a violates 



the convexity of polygon (raqp), while in Figure (b) point q violates the convexity of the 
polygon. 

• Case 2: Both points q and a are outside the triangle A(srp). 

In this case, the edges added to E depend on the angle ^ (srt). There are two cases: 

- Case 2.1: Angle ^(srp) > |. 

In this case, ^(srt) > ^(srp) > |, and therefore, the algorithm adds the edge 
{r, t} to E. Thus, instead of showing 

\sr\ + (1 + V2)\\ra\ + \qp\) + ds^,^,, < (1 + V2)^\spl 

it is enough to show 

\sr\ + (1 + V2)\p\ + 55.,,,, < (1 + V2f\sp\. 

By Corollary [471 



Oq < —\rq\. 

^s,r,q \ ^1 



(7) 



By Observation 



4.2 



(rqp) > and by applying Claim 4.9 on the triangle 



A (rqp) with d-i^we get 



(1 + V2f\qp\ + ^|r(7| < (1 + V2f\ 



Thus, 



\sr\ + (1 + V2f \qp\ + 5s,^,._^ <0 \sr\ + {l + V2f\qp\ + '^\rq\ 



(8) 



\sr\ + 

\sr\ + (l + V2)^{\rr'\ + \pr'\) 
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Inequality (**) follows from triangle inequality for any point r' , and thus it also 



holds for a point r' on {s,p}, such that \sr\ - \sr'\. Therefore, from Claim 4.10 
and since (1 + \/2)^ > ]T:2sin(^i78) ' i^iequality (***) follows. 

Case 2.2: Angle ^{srp) < |. 

Since a and q are outside A(rps), either angle ^{raq) or angle ^{aqp) (facing 
towards s) is less than vr. Assume w.l.o.g. that ^ (raq) < it, thus, by Observa- 
tion 4.2, ^{raq) > Applying Claim 4.9 on triangle ^{raq) with d = 5 gives 



us 



(1 + ^/2)Va| + 2 I"*?! ^ (1 + ^^Z^)^! 



(9) 



Let q' be a point on the intersection of disk Ds^r,p and the extension of {s, q] 
(see. Figure [T]). Then, by convexity we get 



\rq\ + \pq\ < \rq'\ + \pq'\ < 



'I A*) 



\pr\ 



cos(a/2) 



(10) 



The last inequality (*) is obtained by Claim 4.13 



Let 6 be a point on the extension of {s,r}, such that \sb\ - \sp\; therefore. 
By the law of sines, 



\pb\ - \sp\ 



sm{a) I ,sm(a) , . . 

" \^P\ 7^ = 2 sp sm(Q/2). 

cos(f ) 



'sin(f-f) 

Since angle ^(srp) < tt/2, it follows that angle ^{hrp) > 7r/2, thus, 

\pb\ > \pr\. 



(11) 



(12) 
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\sr\ + 5G{r,p) < 



Now we are ready to bound the length of the path: 

sr\ + (1 + \/2)^{\ra\ + \pq\) + '^\qo-\ 
\sr\ + (1 + V2)'^(\rq\ + \pq\) 

cos(a/2) 

cos(a/2) 

cos(a/2) 
sp\{l + (1 + \/2)2(2tan(a/2))) 
sp\{l + (1 + \/2)2(2tan(7r/8))) 
sp\{l + 2(1 + V2){1 + V2){V2- 1)) 
sp|(l + 2(l + \/2)) 
sp\{l + V2f. 

The last inequality (**) follows from the fact that tangent is a monotone in- 
creasing function in the range (0,7r/4]. 




Figure 7: Illustrating the proof of Lemma 4.11 case 2, when ^{srp) < f 



Claim 4.13. Let a, b, and c he three points on a circle, such that ^{abc) - tt - a. Then, 



\ab\ + \bc\ < 



cos(f)- 



Proof. Let /3i be the angle between ba and ca, and let (32 be the angle between be and ca, 
as depicted in Figure |8| By the law of Sines we have 

\bc\ \ab\ 



ac 



sin(7r-a) sin(/3i) sin(/32) 



Therefore, 



I ,i |, I I I sin(/32) sin(/3i) 
|a6| + |6c| = \ac\—, r + \ac\ 



\ac\ 



sin(7r - a) sin(7r - a) sin(7r-a)) 



(sin(/32) + sin(/3i)). 
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For < a < ir/2, this function is maximized when /3i = /32 = ^. Thus, 

sin(f) 



\ab\ + \bc\ < 2\ac\ 
= 2\ac\ 



sin(7r - a) 
, sin(a) 



2\ac\ 



2cos(f) 



ac 



cos(f) 




Figure 8: Ihustrating the proof of Claim 4.13 



Theorem 4.14. For every set of points P, there is a strong planar t-spanner with t = 
(1 + ^2)^ -(5, where 6 is the stretch factor of Delaunay triangulation with bounded degree 7. 



5 Linear-running-time distributed algorithm for planner t- 
spanner with bounded degree 7 

In this section we show a hnear-running-time distributed algorithm that, given a Delaunay 
triangulation of set of points P and clockwise order of the edges around each point p ^ P, 
computes a planner i-spanner of P with bounded degree 7, where the stretch factor t is as 
before {{1 + V2f6). 

The distributed algorithm chooses the edges according to the same principle as in 
Algorithm [T| meaning Lemma [4. 8| still holds, only this time the edges that could be chosen 
by the algorithm are predefined. In order to predefine the candidate edges, we need to 
observe which kind of edges in a cone could be chosen by the algorithm. Moreover, we 
need to compute and store these edges in a sorted order by length in linear time. The 
linear running time is achieved by selecting the edges based on several properties rather 
than sorting them. For simplicity of presentation we assumed that all edges in a cone 
are of different length, however even if this is not the case, we can solve it using known 
method in distribute computing such as wait and notify. Notice that Step |6] in Algorithm [5] 
terminates, since we are in the distributed setting. Moreover, each time-step there exist at 
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least two points that remove at least one edge from the top of their lists. Since there are at 
most linear number of edges that can be removed, and each edge removal takes constant 
time, the total running time of the algorithm is 0{n). 

Algorithm 3 BoundSpanner(P) 

Input: DT(P) - Delaunay triangulation of set of points P and a clockwise order of the 

edges around each point 
Output: A planar f-spanner G - (P,E') with maximum degree 7 

1: Initialize Cp for each p e P /* with respect to edge {p,qmin) */ 

2: for every cone Cp^ in Cp do 

3: CandidateEdgesInCone(Cpj.) /* calling a subroutine that returns a set of edges 

in a cone that are candidates to be added to the 
set E' in a sorted order by length */ 

4: E' ^ 

5: E* ^0 /* additional edges to be included in E' in the end */ 

6: Let List{p) be the sorted list obtained by merging the lists List{Cpj) for all Cp^ e Cp 
7: for each p in P (* in distributed behavior *) do 
8: DistributedEdgeSelction(p) 
9: E' <-E'u E* 



Lemma 5.1. Let {s,pi}, {s,pfc} all the edges of a cone c^Cp in clockwise order. Let 
W*p_^ p_ be the maximum wedge such that \spm-i\ ^ \spm\ for every 1 < m< i. Symmetrically, 
let W*^p.^p^ be the maximum wedge such that \spm+i \ < \spm\ for every j <m < k. I.e., wedge 
W*s,pi,pi (alternatively, W*s^pj,p^) is the largest increasing sequence of edges clockwise 
(counterclockwise) starting from {s,pi} ({s,pk}, respectively). Let Eshort be the set of 
shortest edges in the wedge Wg^p^^pj, i.e., Vei e Escort o^n-d Ve2 e Wg^p^^p^ it holds that 
< |e2|. Then, the edge of cone c that has been chosen by Algorithm^ to be added to E' , 
is from the set E short u ^/,pi,p, ^ 

Proof. Assume on the contrary that Algorithm [s] chooses an edge {s,pt} e c to be added 
to E' that is not in u W^*p^ p. ^W*.p.p^. W.l.o.g. (due to symmetry) assume {s,pt} e 

Ws,p„pr\{{siPr]] for {s,pr} e Eshort- Observe {s,pj}, due to the maximality of W^p.^p^, we 
have \spr\ < \spj\. Since {s,pr} is the shortest edge in the wedge Ws^p^^p^ and since {s,pt} i 
Eshort, it implies that \spr\ < \spt\. Therefore, \spr\ < min{\spt\,\spj\} in contradiction to 



Lemma 4.8 (see, Figure^. 



Lemma 5.2. The resulting graph G - {P,E) of Algorithm^ and the resulting graph G' - 
{P,E') of Algorithmic are identical. 

Proof. Note that the initialization of the cones for each point is the same in both algorithms. 
Moreover, similarly to Algorithm [T| Algorithm [3] chooses an edge from every cone to be 
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Algorithm 4 CandidateEdgesInCone(cone c) 
Input: A cone c in Cp 

Output: A sorted list List{c) by length of edges in the cone c that are candidates to be 
added to the set E' 

1: Let {{p,Pi}, ■ ■ ■ ,{p,Pk}} be the edges in the cone c in clockwise order 

2: Initialize Listi(c) *- {{s,pi}} and List2{c) {{s,pk}} 

3: 

4: while \{p,Pi}\ > \{p,Pi-i}\ do 
5: Add {p,Pi} to Listi{c) 
6: i <- i + 1 
7: i k-1. 

8: while \{p,Pi)\ > \{p,Pi+i]\ do 
9: Add {s,pj} to List2{c) 
10: i 

11: List{c) merge(Lzsti(c), List2{c)) /* merging the lists by length */ 

12: Add the shortest edge in c excluded edges in List{c) in the appropriate place in the 
sorted list List{c) /* in case there are more than one edge - add them all */ 



Algorithm 5 DistributedEdgeSelction(p) 
Input: Sorted list List{p) 

Output: The spanner edges in p contributed by p 
1: while {List{p) ^ 0) do 

2: {p,q} <- TOP{List{p)) /* TOP{List{p)) is the first edge in List{p) */ 

3: Remove {p, q] from List(p) 

4: if not VCp. contain {p,q}, Cp. n E' - then 

5: Go back to[T] 

6: while {\TOP{List(q))\ < \{p,q}\) do (nothing) /* distributed behavior ensures 

while-loop termination */ 

7: if (yCq. contain {p,q}, Cg. n E' ^ 0) then 
8: E'^E'u{{p,q}} 
9: Wedge(p, q) 
10: Wedge(g,p) 
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s 



Figure 9: Illustrating the proof of Lemma 5.1 The edges with no potential to be chosen 
by Algorithm [T] are dashed. 



added to E' . By Lemma 5.1 , for every cone Algorithm [s] chooses an edge among all possible 



edges Algorithm [T] could choose for that cone. Not only that, but the edges are chosen by 
the same principle - the shortest edge, whose cones from both sides "agree" on, is added 
to E. Therefore, E^E' and G^G'. ■ 

Theorem 5.3. Algorithm^ has the following properties: 

1. The output graph of the algorithm has a bounded degree 7. 

2. The stretch factor of the output graph is (1 + \/2)'^6, where 5 is the stretch factor of 
Delaunay triangulation. 



3. The running time of the algorithm is linear. 



Proof. 1. The degree bound follows exactly the degree bound of the previous algorithm 
(Algorithm [T]) and thus is at most 7. 



2. By Lemma 5.2, the output graph of Algorithms [T] and [3] are identical. Therefore, 



the stretch factor of Algorithm [3] is the same as the one of Algorithms [T} which is 
{l + ^/2f5. 

3. Steps [l] to [3] in Algorithm [s] require traversal over all points and edges of DT(P). 
Since DT(P) contains 0{n) edges we get 0{n) running time. Merging two sorted lists 
and inserting edges of the same length cost 0{n) time. The distributed subroutine 
requires every point to remove all edges incident to it from their list and wait, in the 
worst case, until all edges have been removed from their lists. Since every edge is 
incident to two points, the overall running time of this subroutine for each point is 
as the number of edges which is 0{n). Since all points perform this subroutine in a 
distributed manner, we get 0(n) running time. 
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